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ABSTRACT 

We selected 27244 nearby, red, giant early-type galaxies (RGEs) from the Sloan Digital Sky Survey (SDSS). In a large 
fraction (53 %) of their spectra the [O III]A5007 emission line is detected, with an equivalent width (EW) distribution 
strongly clustered around ~0.75 A. The vast majority of those RGEs for which it is possible to derive emission line 
ratios (amounting to about half of the sample) show values characteristic of LINERs. 

The close connection between emission lines and stellar continuum points to stellar processes as the most likely source 
of the bulk of the ionizing photons in RGEs, rather than active nuclei. In particular, the observed EW and optical 
line ratios are consistent with the predictions of models in which the photoionization comes from to hot evolved stars. 
Shocks driven by supernovae or stellar ejecta might also contribute to the ionization budget. 

A minority, ~4%, of the galaxies show emission lines with an equivalent that is width a factor of > 2 greater than 
the sample median. Only among them are Seyfert-like spectra found. Furthermore, 40% of this subgroup have a radio 
counterpart, compared to ~6% of the rest of the sample. These characteristics argue in favor of an AGN origin for their 
emission lines. 

Emission lines diagnostic diagrams do not reveal a distinction between the AGN subset and the other members of the 
sample, and consequently they are not a useful tool for establishing the dominant source of the ionizing photons, which 
is better predicted by the EW of the emission lines. 
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1. Introduction 

The presence of ionized gas in early-type galaxies, part of 
a complex multiphase interstellar medium (ISM), has now 
been firmly established. Among other issues, the source of 
ionization that powers the observed emission lines is still a 
matter of debate. Several alternatives have been advanced, 
ascribing it to active galacti c nuclei (AGN) to hot evolved 
stars , and to shocks (e.g. HeckmanI Il980l : iBinette et al.l 
119941: iDopita fc Sutherlandlll995D . A better understanding 
of this phenomenon is important for several reasons, e.g. 
in the contexts of the origin and evolution of early-type 
galaxies and of their ISM and of a correct census and de- 
mography of active galaxies. 

In this paper we explore this issue by study- 
ing homogeneous sample of galaxies, restricted to a 
rather narrow range of galaxy properties. We isolated 
nearby, quiescent (from the point of view of star for- 
mation) , inassive early-type galaxies from the SDSS 
(jYork et al.l 12000*) . More specifically, from the main sam- 
ple of ~800000 galaxies with spectra available from 
the SDSS, Data Release (DR) 7, we used the MPA- 
JHU DR7 release of spectrum measurements, available 
at http://www.mpa-garching.mpg.de/SDSS/DR7/, to se- 
lect all sources with the following spectroscopic criteria: 
1) redshift between 0.01 and 0.1, 2) Ca break strength 
Z3„(4000) > 1.7, and 3) stellar velocity dispersion cr* > 156 
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km s . If adopting the scaling relation between mass and 
velocity disp ersion this corresponds to galaxies with > 
5 X IO^^Mq, (jHvde fc Bernardil2009fl . We only retain early- 
type galaxies, i.e. objects with a concentration index (the 
ratio of the radii including 90% and 50% of the r bari d 
light) Cr > 2.86 (jNakamura et al.ll2003l: IShen et al.ll2003D . 
This selection yields 27244 red, giant, early-type galaxies 
(hereafter RGEs). 



2. Emission lines and stellar continuum 

We considered the [O III] A5007 emission line flu x measured 
after subtracting a starlight template fsee Kauffmann et al.l 
l2003l for a detailed description of the method used for the 
continuum subtraction). The [O III] line is detected at a 
significance higher than 3a in ~ 53% of RGEs of the sam- 
ple. Figurc[T] compares the [O III] flux with the k-corrected 
i band magnitude within the SDSS fiber, with both quanti- 
ties corrected for galactic absorption. A strong connection 
emerges, and the vast majority of the objects are clustered 
in a narrow str i pe. W e use d the algorithms propos ed by 
iDempster et all (|1977^ and iBucklev fc James! (|l979fl . im- 
plemented as emmethod and buckleyjames in IRAF, to deal 
with censored data in order to derive the best-fitting loga- 
rithmic slope of the F^q mj versus i relation. The resulting 
values are 0.41 ± 0.01 and 0.40 ± 0.01, respectively, both 
consistent with a constant ratio between emission line and 
stellar continuum. 
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Fig. 1. Left panel: logarithm of the [O III] emission line flux (in units of 10~^^ergs^^ cm~^) versus the k-corrected i band 
magnitude within the SDSS fiber, both quantities corrected for galactic absorption. Green dots mark upper limits in line 
flux. The dashed red line corresponds to a constant ratio between the two quantities. Contours (blue for the objects with 
an [O III] detection, yellow for upper limits) represent the iso-densities of DR7 galaxies in the same redshift range; levels 
are in geometric sequence with a common ratio of 4. Right panel: histogram of the residuals from the median line. The 
dashed histogram is the contribution of upper limits. The dispersion of the distribution is 0.18 dex. Bottom: residuals 
from the median vs redshift. 



In the right panel we show the histogram of the resid- 
uals from the median line that is reproduced well by a 
Gaussian distribution with a dispersion of only 0.18 dex 
(similar results, although with slightly broader distribu- 
tions, are obtained using different lines and/or continuum 
bands). Nonetheless, there are outliers with a line excess, 
R[0IIlI3, of up to 100. Setting a threshold at R[OIII] = 
5, we find 0.5% outliers, a fraction that increases to 1.2% 
(3.5%) when lowering the limit to 3 (2). 

The EW of the [O III] line, estimated by considering the 
continuum level 200 pixels around the line, has a median 
value of ~ 0.75 A including upper limits in the analysis) 
and the EW dispersion is of 0.20 dex; 95 % of RGEs with 
an [OIII] detection have 0.3 < EW < 1.7 A. 

Contours in Figur£[T] represent the iso-densities of the 
~340000 DR7 galaxies in the same redshift range, 0.01 < 
z < 0.1, but without any further selection of the spectro- 



^ Defined as the ratio between the [O III] flux of a galaxy with 
respect to its median value at a given magnitude. 



photometric properties. The general galaxies population 
shows a concentration in the same locus as covered by the 
RGEs, but also a substantial population extending toward 
high line excesses (^ 30% of the sample has an excess of 
R[0III]>3). 

We looked for trends between the residuals and the spec- 
trophotometric parameters used for the sample selection 
(namely redshift, I?„(4000), concentration index, and ve- 
locity dispersion) , but we failed to find any statistically sig- 
nificant link. In particular, there are no differences between 
RGEs at low and high redshifts: the best linear fit of the 
residuals against redshift indicates that objects at z=0.01 
have a ratio between line and continuum only (20 ± 11)% 
higher than those at z=0.1. 



3. Spectroscopic diagnostic diagrams 

Fig. [2] shows the l ocation in the spectros copic diagnos- 
tic diagrams (e.g. HeckmanI 1980: Baldw in et al.l [l981l : 
IVeilleux fc Osterbrockl Il987t Ikewlev et al.i 120061 ) for the 
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Fig. 2. Spectroscopic diagnostic diagrams for RGEs. The solid lines are from iKewlev et al] (|2006f) and separate star- 
forming galaxies, LINER, and Seyfert; in the first panel the region between the two curves is populated by the composite 
galaxies. Blue dots mark the "strong" outliers (objects with a line excess with respect to the median value R[0III]>5) 
in Figure[Tl and red dots the "weak" outliers (3<R[0 III]<5). Contours represent the iso-densities of all DR7 emission 
line galaxies. 
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Fig. 3. [N II]/Ha ratio for all sources with both lines de- 
tected. Emission lines of the objects above the horizontal 
dashed line are unlikely to be powered by star formation. 



RGEs that have all the relevant emission lines detected at 
SNR>3 separately for each diagram. Starting from the left 
side, the percentages with respect to the whole sample are 
of 29, 23, and 8% in the three diagrams, respectively. The 
vast majority of the objects fall in the LINERs region, while 
the Seyfert and star-forming regions are scarcely populated. 

Considering the properties of outliers first, 'strong' out- 
liers (i.e. the objects with R[0 III] > 5) represent essentially 
all of the Seyferts; nonetheless, most of them are LINERs. 
'Weak' outliers (with 3 < R[OIII]< 5) do not differ sig- 
nificantly, on the one hand, from the 'stronger' ones (but 



there are fewer Seyfert among them), but on the other, they 
cannot be readily separated from the bulk of the RGEs pop- 
ulation. 

We now turn our attention to the objects that have no 
optical classification because at least one of the key emis- 
sion lines is undetected. This subsample is clearly very im- 
portant, because it is composed of ~ 3/4 of the galaxies. 
Nonetheless, [N II] and Ha are both detected in ^50% of 
the galaxies and this allows us to derive at least a tenta- 
tive classification. In fact, log([N II]/Ha) > -0.3 in 90 % 
of them (see Figur£|31), a threshold above which no star- 
forming galaxies are found. On the other hand, Seyfert 
galaxies usually have bright emission lines, thence they are 
expected to be all properly cataloged by the diagnostic di- 
agrams in Fig. [5] This leads to the conclusion that these 
galaxies, amounting to about half of the RGEs, are gener- 
ally LINERs. 

4. Error budget of the emission line fluxes 

Since we are studying emission lines generally of very low 
EW, it is important to assess the reliability of such mea- 
surements and of the error estimates. Two sources of uncer- 
tainties should be considered, the first related to statistical 
errors, the second to the effects of starlight modeling. 

We deal with statistical uncertainties by studying the 
galaxies with multiple SDSS spectroscopic observations. 
There are 3103 such objects: 45% of them have the [O III] 
line detected in both datasets, 35% are upper limits in 
both observations, while 20% change from detected to un- 
detected (or vice- versa). Considering the objects detected 
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Fig. 4. Difference between pairs of duplicated [O III] flux 
measurements divided by the error estimated propagating 
the individual uncertainties. The solid line is a Gaussian 
distribution with a dispersion of 1. 



in both observations, the differences between the indepen- 
dent measurements closely follow a Gaussian distribution 
with a dispersion equal to what is expected based on the 
uncertainties in each flux estimate (see Fig S]). We con- 
clude that the error assessment is robust and that most 
(85%) of the [O III] detections are confirmed by multiple 
observations. 

In addition to the statistical uncertainties, we must also 
consider the effects related to the accuracy of the sub- 
traction of the stella r emission. As discussed in detail by 
lAnnibali et al.l (|2010t ). because of the degeneracy between 
age, metallicity, and extinction, fits to the continuum emis- 
sion of similarly good quality can produce significantly dif- 
ferent residual spectra. This is particular relevant for the 
H/3 line (and to a lesser extent for Ha) which is superim- 
posed on the stellar B almer a bsorption feature. The method 
proposed by Annibal i et al.l ([2010) to cope with this prob- 
lem in their study of a sample of nearby ETG is to explore 
the changes in the emission lines fluxes related to the var- 
ious stellar templates, limiting it to those that yield a fit 
with < 2. The adopted emission line flux is a weighted 
average of those derived from all acceptable fits, while its 
error is a weighted value of the dispersion between the in- 
dividual estimates. 

The very large number of SDSS objects considered pre- 
vents a detailed analysis of the error budget considered by 
these authors. However, this provides us with guidance on 
the typical uncertainties that should be associated with the 
line measurements due to a template mismatch. From their 
published values for their smaller synthetic aperture (com- 
parable in size with the SDSS fiber), we estimated that the 
median errors in the EW for H/3, [O HI], and Ha are 0.16, 
0.08, and 0.08 A, respectively. We adopt these reference 
values also for the SDSS sources. The errors in EW have 
been been converted into a line flux error and then added 
quadratically to the uncertainties provided by the SDSS 
database. 

Our results for the connection between [O HI] flux and 
z magnitude are only marginally changed by adopting this 
more conservative error treatment. The fraction of detected 



sources decreases from 53 to 50%, but the slope of the re- 
lation remains unchanged. The number of outliers is not 
affected. 

We then reconsidered the diagnostic diagrams, exclud- 
ing all sources that do not meet the 3(t criterion with the 
revised errors. The fraction of sources in the various dia- 
grams decreases by 25 to 40 %, but their visual appear- 
ance is essentially preserved. In order to quantify more 
subtle differences, we estimated the changes of median val- 
ues and dispersions of the various ratios. Not surprisingly 
log [QUI] /H/3 is the more affected ratio, because it changes 
from 0.15 (with an rms of 0.20) to 0.09, (rms = 0.18 dex). 
This is most likely due to the preferred exclusion of galax- 
ies with low H/3 EW. Nonetheless, the result that the vast 
majority of the sources are LINERs is confirmed. 

Finally, we explored how crucial the precise value of the 
adopted errors on EW is. Doubling the reference values, 
the number of sources in the diagnostic diagrams keeps de- 
creasing, but their location is effectively unchanged, with a 
median of log [O III]/H^ = 0.05. 

5. Discussion 

Most nearby, red, giant early-type galaxies are emission line 
galaxies. In particular, in 53% of their SDSS spectra, the 
[O III]A5007 line is detected. The majority of the objects for 
which this analysis was possible show emission line ratios 
characteristic of LINERs. 

The [O III] flux shows a strong correlation with the flux 
measured within the 3" SDSS fiber in the i band. The i 
band magnitude (being less affected by uncertainties in the 
k-correction, absorption, and by differences in the stellar 
population) is a good estimator of the stellar mass within 
the fiber. Thus the ratio between lines and stellar mass is 
essentially constant (showing a dispersion of only 0.18 dex) 
while both quantities vary by a factor of ~ 30. Furthermore, 
there is only a marginal trend for a change of EW with 
redshift: objects at z=0.01 have a ratio between line and 
continuum (20 ± 11)% higher than those at z=0.1 despite 
the change in the size of the region covered by the SDSS 
fiber from ~0.6 to 5.5 kpc. 

These two results are not straightforward for explaining 
if the emission lines are powered by an AGN. In fact, they 
require a fine tuning between the strength of the nuclear 
ionizing field, the spatial distribution of the emission lines, 
and the stellar mass. A general link between AGN activity 
and the toial stellar mass can be envisaged, via the con- 
s tant ratio between the mass of stars and of the black hole 
(■ Marconi fc Huntll2003h . while we find a connection with 
the stellar mass covered by the SDSS fiber. Furthermore, 
also the accretion rate, the second parameter that together 
with the black hole mass sets the radiative output of the 
AGN, should be closely linked with the amount of stars 
within the central 3" of each given galaxy. 

Conversely, these findings strongly suggest that pro- 
cesses related to the stellar population are at the origin 
of the observed emission lines. 



As argued by sever a l authors (e.g. iBinette et al.l 11994 



ly sever al i 

Macchetto et all 119961 : IStasinska et al.l 120081: ISarzi et al.l 
2010D . hot evolved stars, such as post- asymptotic giant 



branch stars (pAGB) and white dwarves (WD), can pro- 
duce a substantial diffuse field of ionizing photons. More 
quantitatively, the observed median EW for Ha in the sam- 
ple considered here is EWhq 0.8 A, well within the range 
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0.6 - 1.7 A predicted bv iBinette et all (|l994f l for a stellar 
population of age ^ 10^^° years, assuming that the cold gas 
intercepts all the ionizing photons. Stasinska ct al. (2008) 
have explored photoionization models in which the Lyman 
continuum radiation is directly estimated from a stellar 
population analysis, and find that, by varying the metallic- 
ity and ionization parameter U, the resulting emission line 
ratios cover the whole region typical of LINERs in the spec- 
troscopic diagnostic diagrams. In particular, the location of 
RGEs is reproduced well by models with 0.5 < Z^,/Zq < 2 
and log U 3.7. 

Shocks can represent an additional so urce of ionization. 
It has been shown by, e.g., lAUen et al.l ((2008) that they 
can produce emission lines with ratios mimicking those of 
active nuclei and of LINERs, in particular. The connection 
between stellar mass and emission line flux suggests that 
such shocks must be driven by the stellar population, as in 
the case of supernovae or of stellar ejecta, rather than by 
bars of triaxial perturbations. 

We conclude that the properties of emission lines in 
RGEs can in general be satisfactorily accounted for by pho- 
toionization related to hot evolved stars or by shocks driven 
by the stellar processes (but below we discuss the few and 
important exceptions of likely AGN dominated objects). 
The narrow spread of observed EW, very difficult to ex- 
plain in an AGN framework, set strong constraints even 
when adopting a "stellar" interpretation. For example, in 
the scenario of lines powered by an evolved stellar popu- 
lation, the apparently stringent requirement that the cold 
gas must intercept essentially all of the ionizing photons, 
becomes probably instrumental in producing a narrow dis- 
tribution of lines EW. In fact this effectively removes the 
fundamental degree of freedom related to the ISM optical 
depth to ionizing photons. The selection of galaxies within 
a rather narrow range of the various spectro-photometric 
parameters space is also a key ingredient in reducing the 
scatter between F[OIII] and i magnitude (see Figured]). 

These results were obtained by focusing on RGEs, but 
the consequences are probably more far reaching. In fact, 
in the F[OIII] vs. i plane a strong concentration is also 
seen for the general population of galaxies, in the same 
region covered by RGEs. This suggests a wider application, 
possibly extending to less massive early-type galaxies and 
to the bulges of spiral galaxies. 

Nonetheless, there are outliers from the general link be- 
tween stellar and line emission. Only among these sources 
we found Seyfert-like spectra. Furthermore, a strong con- 
nection emerges between the excess in line emission and 
the presence of a nuclear radio source, a characteristic sig- 
nature of an AGN. Our selection criteria isolate galaxies 
whose properties are closely matched to the typical hosts 
of low-luminosity radio-loud AGN (Baldi & Capetti 201C|). 
Considering the outliers with R[OIII]> 5 (3), 49% (45%) 
have a radio counterpart in the FIRSlQ catalog within 2" 
from the optical position. In the rest of the sample, only 
6% of the sources have a radio counterpart. The connection 
between outliers and radiosources is preserved even at very 
small line excesses, since 31% of those with 2 <R[0 III]< 3 
and 17% of the objects with an excess by only a factor of 1.5 
to 2, are associated with a FIRST source. Thus, the galax- 
ies that exceed the rather narrow range of line emission EW 



^ Faint Images of the Radio Sky at Twenty centimeters survey 
ijBecker et al.lll995l ) 



set by the stellar processes are very often associated with 
an active nucleus. The observed line excess is likely to stem 
from the additional source of ionizing photons represented 
by the AGN. The rather loose relation between radio power 
and line luminosity seen in these low-luminosity AGN (see 
the results presented bv lBaldi fc Capettil 1201011 might ex- 
plain the incomplete coincidence between line excess and 
the presence of a radio source and, on the other hand, the 
reason not all radio sources produce a detectable increase 
in the line EW. 

Apparently, the location in the spectroscopic diagnos- 
tic diagrams is not sufficient to establish the dominant 
source of ionizing photons, which, conversely, is best pre- 
dicted by the EW of the emission lines. In fact, the likely 
AGN show line ratios that are essentially indistinguishable 
(leaving aside the few Seyfert-like objects) from those prob- 
ably d ominated by ste l lar pr ocesses. As already pointed 
out by IStasinska et al.l (|2008f ) , this has important ramifi- 
cations where the census and the properties of AGN are 
concerned. To isolate and explore the properties of genuine 
low-luminosity active galaxies, spectroscopic observations 
obtained with a smaller aperture are needed to increase the 
contrast between AGN and stellar induced line emission. 
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